The rapidity and centrality dependence of nuclear modification factors at RHIC — 
what does bulk particle production tell us about the nuclear medium? 
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The BRAHMS experiment at RHIC has measured the production of charged hadrons as a function 
of pseudorapidity and transverse momentum in Au + Au, d + Au and p + p collisions at a common 
energy of y / sjvjv = 200 GeV, and from these spectra we construct the nuclear modification factors 
for both "hot" and "cold" nuclear matter. 

In this contribution I will show how these factors evolve with pseudorapidity and collision cen- 
trality. We see a Cronin-like enhancement in d + Au collisions at midrapidity, going to a strong 
suppression at r\ > 2. In central Au + Au collisions we find a suppression both at mid- and for- 
ward rapidities that vanishes for peripheral collisions. We interpret this as signs of several different 
medium related effects modifying bulk particle production in Au + Au and d + Au collisions at RHIC 
energies. 



The processes that govern bulk particle production 
in ultrarelativistic nucleus-nucleus collisions jj are not 
well understood. Exploring this experimentally is one of 
the goals of the RHIC accelerator at Brookhaven Na- 
tional Lab., and of the BRAHMS experiment in par- 
ticular. RHIC has now collided gold nuclei on gold 
(Au + Au), protons on protons (p + p) and deuterons 
on gold (d + Au) at the same center-of-mass energy of 
v/ sjvjv = 200 GeV, allowing us to compare the complex 
dynamics of a heavy-ion collision to the comparatively 
simple case of a nucleon-nucleon interaction, and the 
particle production mechanisms in hot nuclear matter 
(central Au + Au) to those in a cold nuclear medium 
(d + Au). 

The hope is that clear, systematic effects will emerge 
from these comparisons that will tell us about the gross 
properties of the system created in these collisions. Be- 
low we present BRAHMS data on charged hadron pro- 
duction from all three colliding systems, and construct 
the nuclear modification factors RauAu and RdAu that 
are sensitive to such medium-related effects. We will 



show that bulk properties do indeed emerge when com- 
paring both extended systems to p + p collisions, and 
indicate how this can be interpreted in terms of several 
distinct mechanisms for modification of particle produc- 
tion in different parts of r\ — p? space. Here pt is the 
particle transverse momentum, r\ = — ln(tan tj) its pseu- 
dorapidity and 9 the angle of the particles momentum 
relative to the beamline. 

The BRAHMS experimental setup 0, El , one of four 
detectors at RHIC, is a two-arm hadronic magnetic spec- 
trometer with three main parts - see Figure ^ The 
first is a set of event characterization detectors including 
Cherenkov scintillator arrays for multiplicity and colli- 
sion point determination (BB), silicon strips and plas- 
tic scintillator tiles for particle multiplicity and central- 
ity determination (SMA & TMA), and two zero-degree 
calorimeters (ZDC) used for collision point, centrality 
and cross section determination. 

Secondly there is a midrapidity spectrometer arm 
(MRS) that can be positioned at angles from 90 to 30 
degrees with respect to the beam axis. The MRS con- 
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FIG. 1: The BRAHMS experiment at RHIC. For details, refer 
to 0. 



sists of two TPCs (TPM1, TPM2) and a dipole magnet 
for tracking and momentum determination, and a time- 
of-fiight wall (TOFW) for particle identification. The 
MRS can identify tt* and up to p = 2.0 GeV/c and 
p/p up to p = 3.0 GeV/c. At a single angle and magnetic 
field setting the solid angle coverage of the MRS is quite 
limited, but by rotating the spectrometer and changing 
the field of the bending magnet a large coverage in ra- 
pidity and transverse momentum is achieved. 

Finally, the part that makes BRAHMS unique at RHIC 
is the forward rapidity spectrometer arm (FS), which ro- 
tates from 20 degrees (front part only) to 4 degrees rela- 
tive to the beam axis. It consists of two TPCs and three 
Drift Chambers (T1-T5), four dipole magnets (D1-D4), 
one threshold and one ring-imaging Cherenkov counter 
(CI and RICH). The full FS allows identification of pi- 
ons, kaons and protons up to 20 GeV/c. As for the MRS, 
the limited geometrical acceptance is extended by rotat- 
ing the FS and changing the fields. All in all BRAHMS 
can identify pions from midrapidity and up to y ~ 3.5, 
with the RHIC beam rapidity at yt> = 5.3. 

For the p + p and d + Au data sets taken in 2003, the 
above setup was augmented by two sets of plastic scintil- 
lator rings around the beampipe (not shown in Fig.jTj see 
Ref. [J for a description), four rings on each side of the 
nominal interaction point. These counters provided event 
triggering, interaction point and cross section measure- 
ments in the relatively low multiplicity events from these 
collision systems. For further details of the BRAHMS 
experimental setup, see Ref. Q- Also note that for the 
d + Au data presented here, the spectrometers saw only 
the deuteron fragmentation side of the collision and we 
consequently define this to have positive rapidity. 

Figures El and rj 5, 6] show BRAHMS charged hadron 
spectra from Au + Au, d + Au and p + p collisions at 
y/SNN = 200 GeV for a number of pseudorapidities and 
centrality classes. Measuring charged particle production 
with BRAHMS is a multi-step process due to the lim- 
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FIG. 2: Invariant spectra of charged hadrons from Au + Au 
collisions at y^ijvjv = 200 GeV for pseudorapidities r\ = (left 
panel) and r\ = 2.2 (right panel). Various centrality cuts are 
shown for Au+Au. The p + p reference spectra (appropriately 
acceptance scaled) are shown for comparison. The left panel 
also shows the d + Au spectrum. For clarity, some spectra 
have been divided by the indicated factors. Published in [f|. 



ited acceptance — total py spectra of charged hadrons 
at one pseudorapidity |ESl3 as shown are constructed 
by combining several such settings. Our spectra have 
been corrected for finite geometrical acceptance as well 
as trigger and detector efficiency using a GEANT-based 
Monte Carlo simulation of the entire detector setup. De- 
tails of the analysis can be found in Refs. 123 • 

Our aim is to extract any effects of a nuclear medium 
on the total particle production, by comparing to elemen- 
tary nucleon-nucleon collisions where we do not expect 
such a medium to be produced. We construct the Nu- 
clear Modification Factor Rab- 



Rab 



1 N AB {p Tl ri) 



(1) 



where A and B are the colliding systems (i.e. Au or 
d) and (N^B) i s the average number of binary nucleon- 
nucleon collisions in such an event as estimated by the 
Glauber model. Nab and N pp are the charged hadron 
productions in the A+B and p + p collisions respectively, 
taken at the same transverse momentum px and pseudo- 
rapidity r\. 

If there were no medium effects, i.e. if a nucleus- 
nucleus collision could be viewed simply as individual 
nucleon-nucleon collisions, then we would expect an Rab 
factor of unity above a certain pt threshold and decreas- 
ing smoothly below. The latter is because at low px we 
expect a scaling of particle production with the num- 
ber of participating nucleons (N part } rather than (N co u). 
Both for Au + Au and d + Au collisions, we see strong 
deviations from this simple behavior. 
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FIG. 3: Invariant yield distributions for charged hadrons produced in d + Au and p + p collisions at -*/s = 200 GeV at 
pseudorapidities 7/ = 0,1.0,2.2,3.2. Published in Q. 
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FIG. 4: Top row: Nuclear modification factors RauAu as a 
function of transverse momentum for Au + Au collisions at 
r\ = and r\ = 2.2 for the — 10% most central collisions. 
Middle row: as top row, but for centralities 40 — 60%. Bot- 
tom row: ratio of the RauAu factors for the most central and 
most peripheral collisions at the two rapidities. The dotted 
and dashed lines show the expected value of RauA u using a 
scaling by the number of participants and by the number of 
binary collisions, respectively. Error bars are statistical. The 
gray bands indicate the estimated systematic errors. The gray 
band at pr — is the uncertainty on the scale. Published 
in g. 



In central Au + Au collisions, where we expect to pro- 
duce "hot" nuclear matter (see e.g. @), we see a strong 
suppression of particle production relative to p + p scal- 
ing — see Fig. H 5], upper left panel. The expected rise 
at low pt is seen up to pr ~ 2 GeV/c, but then RauAu 
seems to level off at a value well below unity and subse- 



quently decrease again. This behavior is in sharp contrast 
to what was seen at SPS energies, i.e. Pb+Pb collisions at 
center-of-mass energies of up to ^snn = 17 GeV, where 
one observed a strong enhancement of Rptpt, dubbed 
the Cronin enhancement |l0j| . interpreted as an initial 
state broadening of the distribution of quark momenta 
(k,T broadening). We also see (middle left panel of Fig.^J 
that for more peripheral Au 4- Au collisions the suppres- 
sion vanishes and we reach the behavior expected from 
p + p scaling. This indicates that we may have an on- 
set of a new mechanism for suppressing high-py parti- 
cle production between top SPS and RHIC energies, and 
between peripheral and central collisions at RHIC — i.e. 
in the collisions where we reach the highest energy densi- 
ties |23j. The behavior reported here has also been seen 
by the other RHIC experiments [Tl|. 

To extend this picture, BRAHMS has also measured 
RauAu from central collisions at a more forward rapid- 
ity i] w 2.2, shown in the upper right panel of Fig. 0] 
The suppression seen at midrapidity is still present here, 
showing that if it is due to some final-state mechanism 
prevalent in the fireball, the medium that has this sup- 
pressing property must extend out at least two units of 
pseudorapidity away from rj = 0. Again, the suppression 
vanishes for more peripheral collisions. 

Using the above definition of the nuclear modifica- 
tion factor we need to compare separate measurements 
of nucleus-nucleus and p + p collisions, introducing ad- 
ditional systematic errors. To avoid this problem, we 
can approximate p + p collisions by peripheral nucleus- 
nucleus collisions scaled to ^/{N^ l " ph ) and rather con- 
struct the ratio of central to peripheral collisions Rcp- 
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The bottom panels of Fig. 0] show the ratio Rcp for 
Au + Au collisions, using the centrality classes 0-10% and 
40-60%. Again we see a clear saturation of the modifi- 
cation factor at a value below unity and a subsequent 
further drop as pt increases, showing that there is an 
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FIG. 5: Nuclear modification factor for charged hadrons at pseudorapidities n = 0,1.0,2.2,3.2. One standard deviation 
statistical errors are shown with error bars. Systematic errors are shown with shaded boxes. The vertical shaded band around 
unity indicates the error on the normalization to (Ncoll)- Published in |(|. 



attenuating mechanism present in the central collisions 
that is either not there or greatly reduced in peripheral 
collisions. 

Figure shows the corresponding ratios RdAu for 
minimum-bias d + Au collisions at v / s]v7V = 200 GeV, for 
pseudorapidities r\ = 0, 1.0 , 2.2 and 3.2 [(|. From these 
collisions at midrapidity, where the size of the fireball is 
smaller and we don't expect to produce an extended hot 
and dense medium, we see an enhancement of RdAu very 
similar to the Cronin effect, rather than a suppression. 
However as we move to more forward rapidities, this ef- 
fect vanishes and at r\ > 2.2 we once again see indications 
of a suppression. 

In Figure [5] (j we have defined three centrality classes 
(0-20%, 30-50% and 60-80%) and again constructed 
the ratio Rcp of central (or semi-central) to peripheral 
collisions. While the factors are approximately equal 
at midrapidity, there is a clear change of behavior to- 
ward higher rapidities as the Rcp for central collisions 
drops significantly below that for semi-peripheral colli- 
sions. This indicates that we have an attenuating effect 
on forward particle production when the deuteron hits 
the Lorentz contracted gold nucleus well away from the 
edges, i.e. where we have a constant, dense nuclear initial 
state. 

To interpret these observations, several distinct mech- 
anisms may be needed. Firstly, the midrapidity enhance- 
ment of RdAu can be explained through a broadening of 
the initial momentum distributions due to multiple scat- 
tering (Cronin enhancement) as for the SPS data, indi- 
cating similar initial state medium effects at midrapidity 
in d+Au collisions at ^/snn — 200 GeV and in in Pb+Pb 
collisions at ^snn — 17 GeV. 

Secondly and more interestingly, the suppression of 
high-pT particle production in central Au + Au colli- 
sions is indicative of partonic energy loss due to gluonic 
bremsstrahlung in a colored medium. In other words, 
the partons lose energy through gluonic interactions be- 
cause they traverse a medium with a high density of color 
charges. The fact that the suppression disappears in 



more peripheral collisions and turns to a Cronin-like en- 
hancement for central or min-bias d+Au collisions where 
we do not expect to produce a deconfined medium, indi- 
cates that this is indeed a final-state effect and not just a 
consequence of the high-energy nuclear initial state. Re- 
sults from the STAR experiment on strong suppression 
of far-side jets in central Au + Au collisions extend and 
support this conclusion 

Thirdly, the transition from a Cronin-like enhance- 
ment at midrapidity to a clear suppression with increas- 
ing rapidity in d + Au collisions indicates yet another 
mechanism, related to the nuclear initial state. Recent 
theoretical work has predicted and effect similar to this 
on the basis of gluon-saturation phenomena, a scenario 
known as the Color Glass Condensate [la . fLU fl5j | . It is 
based on the observation from deep-inelastic scattering 
at HERA 

that as the momentum transfer between 
the electron and the struck parton increases, the observed 
gluon density function seems to diverge at small values 
of x, i.e. there is a high density of so-called "wee glu- 
ons" carrying a very small fraction x of the total nucleon 
momentum. Due to the finite size of a nucleon, one can 
assume that at some density these small-a; gluons will 
start to fuse due to gluonic self-interaction. For a nuclear 
system at high energies, one can then predict a transition 
to a condensed state of "colored glass" — a state with 
a high density of color charges that evolves slowly com- 
pared to the timescale of an ultrarelativistic heavy-ion 
collision (« 1/m/c), and with the gluonic wavefunctions 
extended in the transverse directions. The transverse 
momentum scale Q 2 S for the onset of this saturation in 
a nuclear collision can be shown to depend on the gluon 
density and thus on the number of participating nucleons, 
and is connected with the rapidity of measured particles 
by Q 2 S ~ Ah Xy (see e.g. (Tj), where A - 0.2 - 0.3 is 
obtained from fits to HERA data. At RHIC and for the 
Pt range covered by BRAHMS, the onset of saturation 
is predicted [TBI to occur between rapidities and 3, in 
agreement with our experimental observations. 

However, in this part of phase space nuclear shadow- 
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FIG. 6: Central (full points) and semi-central (open points) Rcp ratios from d + Au collisions (see text for details) at 
pseudorapidities 77 = 0,1.0,2.2,3.2. Published in [g. 



ing [l| is also important. The shado wing effect, observed 
e.g. in muon-nucleus interactions [l7j . is a depletion 
of low-a; partons in nuclear systems. It has been in- 
terpreted both as destructive interference reducing the 
probability of interactions with partons at the back of 
an incoming nucleus, and as gluon recombination at high 
densities but without the condensation mechanism dis- 
cussed above. This effect is certainly present in our data, 
as shown by the observation that simulations using the 
HIJING event generator which includes shadowing well 
reproduces BRAHMS data on total chargedparticle pro- 
duction at all accessible pseudorapidities 4]. A recent 
discussion in Ref. 0] indicates that such effects may not 
be enough to describe the results presented here, but this 
discussion relies on extrapolated fits to RHIC midrapid- 
ity data. More data at forward rapidities at higher pt 
and subsequent theoretical work is needed to determine 
whether the observed suppression is indeed partly due to 
the formation of a Color Glass Condensate in the initial 
state at RHIC. Also, a better treatment of soft beam frag- 
mentation effects may be needed to properly construct 
the RdAu factor. 

An effect that may affect the data and that has not 
been discussed in this contribution is a possible isospin 
dependent contribution to the nuclear modification fac- 
tors. A nucleus has both protons and neutrons, and so 
the denominator in the Rab factor should be a properly 
weighted average of p + p, p + n and n + n collisions. 
Data on this at SPS energies have been taken by the 
NA49 collaboration, but are unfortunately not yet pub- 
lished. However, while this may have an effect on the ex- 
act shapes of the Rab factors, the consequences should 
not be large enough to affect the overall trends presented 
above. 

In conclusion, BRAHMS has measured the charged 
hadron production from Au + Au, d + Au and p + p col- 
lisions at the same center-of-mass energy of ^snn = 
200 GeV, both as a function of pseudorapidity and cen- 
trality. Our measurements cover an extensive part of 



phase space, uniquely including the forward rapidity re- 
gion which probes the small-x part of the nuclear struc- 
ture functions. Through the nuclear modification fac- 
tors RauAu and RdAu we have uncovered a suppression 
of high-pT particle production in d + Au collisions at 
forward rapidities, in contrast to the enhancement pre- 
viously seen at midrapidity. We also find a suppression 
of RauAu at high-pr in central collisions, extending at 
least two units of pseudorapidity away from r\ = 0. The 
results presented here can be interpreted in terms of a 
number of medium-related effects modifying bulk parti- 
cle production: 

• A Cronin-like enhancement in d + Au collisions at 
midrapidity, due to initial state broadening of in- 
trinsic quark transverse momenta. 

• A final-state suppression of high-py particle pro- 
duction in central Au+Au collisions, interpreted as 
a partonic energy loss due to gluon bremsstrahlung 
in a deconfincd medium with a high density of color 
charges. 

• An initial state suppression of RdAu at forward ra- 
pidities, at least partly due to nuclear shadowing 
effects depleting the number of small-cc partons 
available for interaction. There may also be an ad- 
ditional suppressing effect from the saturation and 
condensation of small-x gluons, as predicted in the 
Color Glass Condensate scenario. 
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